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Abstract 

We present an updated next-to-leading order analysis of the B — > Xgj branching 
ratio and photon spectrum, including consistently the effects of Fermi motion in 
the heavy-quark expansion. For the Standard Model, we obtain B(i3 ^s7) = 
(2.57 lb 0.26^Q3g) X 10""^ for the integral over the high-energy part of the photon 
spectrum with El^^ > 2.2 GeV, where the first error reflects the uncertainty in the 
input parameters, and the second one the uncertainty in the calculation of Fermi 
motion. This prediction agrees with the CLEO measurement of the same quantity 
within one standard deviation. From a reanalysis of the CLEO data, we obtain 
for the total branching ratio B{B Xs^f) = (2.62 ± 0.60exp toJoth) ^ 10"^ using 
the measured rate above 2.2 GeV, and (2.66 ± 0.56exp 48 th) ^ 10~^ using a fit to 
the photon energy spectrum. Both values are consistent with the Standard Model 
prediction of (3.29 ib 0.33) x 10^^. Our analysis contains an improved discussion 
of renormalization scale dependence and QED corrections. We also study the 
sensitivity of the branching ratio and photon spectrum to hadronic parameters 
such as the 6-quark mass, and to possible contributions from New Physics beyond 
the Standard Model. 
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1 Introduction 



The inclusive radiative decays B Xg'y have been the subject of a considerable number 
of experimental and theoretical investigations. About three years ago, the CLEO Col- 
laboration reported a first measurement of the branching ratio for these decays, yielding 

i 

B{B -> X,7) = (2.32 ± 0.57 ± 0.35) x 10"^ , (1) 

where the first error is statistical and the second is systematic (including model de- 
pendence). Recently, the ALEPH Collaboration has reported a measurement of the 
corresponding branching ratio for 6-hadrons produced at the Z resonance, yielding 

B{Hb X,7) = (3.11 ± 0.80 ± 0.72) x 10"^ , (2) 

which is compatible with the CLEO result. 

Being rare processes mediated by loop diagrams, inclusive radiative decays are po- 
tentially sensitive probes of New Physics beyond the Standard Model, provided a precise 
theoretical calculation of the branching ratio can be performed. The general framework 
for such a calculation is provided by the heavy-quark expansion, which predicts that, 
up to small bound-state corrections, inclusive decay rates agree with the parton model 
rates for the underlying decays of the b quark Q-f^. As long as the fine structure of 
the photon energy spectrum is not probed locally, the theoretical analysis of -B ^ Xsj 
decays relies only on the weak assumption of global quark-hadron duality. The leading 
nonperturbative corrections have been studied in detail and are well understood 
Still, the theoretical prediction for the branching ratio suffers from large perturbative 
uncertainties of about 30% if only leading-order expressions for the Wilson coefficient 
functions in the effective weak Hamiltonian are employed [0-[|T^. Therefore, it was 



an important achievement when last year the full next-to-leading order calculation of 
the total B —>■ Xgj branching ratio in the Standard Model was completed, combining 
consistently results for the matching conditions ||15| , |16|, matrix elements |T^, [18| and 



anomalous dimensions [O. With this calculation the theoretical uncertainty was re- 



duced to a level of about 10% [T^, More recently, the next-to- leading order analysis 
was also extended to two-Higgs-doublet models ||2l|, 

Whereas considerable effort has thus gone into calculating the total B Xs'~f branch- 
ing ratio, little progress has been made in understanding the structure of the photon 
energy spectrum at next-to-leading order. On the other hand, what is experimentally 
accessible is only the high-energy part of the photon spectrum, and an understanding 
of the spectral shape is thus a prerequisite for extrapolating the data to the full phase 
space. For instance, the CLEO Collaboration has measured the spectrum in the energy 
range between 2.2 and 2.7 GeV (in the laboratory) and applied a correction factor of 
0.87 ± 0.06 to extrapolate to the total decay rate p3[.[] This factor does not take into 
account the full next-to-leading order corrections to the decay rate. More importantly. 



it relies on an estimate of bound-state effects |Q obtained using the phenomenological 



similar treatment is followed in the ALEPH analysis 
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model of Altarelli et al. , which is not fully consistent with the heavy-quark expansion. 
The small uncertainty assigned to the correction factor reflects the fact that the model 
parameters (the Fermi momentum and the constituent quark masses) have been tuned 
to fit the lepton spectrum in 5 — > Xc £ z/ decays and then used to predict the photon 
spectrum in B Xg'j decays. It is now known that there is no theoretical justification 
for such a treatment |2^-[^. These conceptual shortcomings were not improved in the 
updated analysis of the photon spectrum presented by Ali and Greub a few years ago 
ir?! , although more complete formulae for the perturbative corrections were used in this 
work and a more conservative error analysis was presented. 

The fact that only the high-energy part of the photon spectrum in B ^ Xs'~f decays 
is accessible experimentally introduces a significant additional theoretical uncertainty 



29| , which has been ignored in previous analyses. This observation limits the potential 



of existing data on these decays to probe or constrain New Physics beyond the Standard 
Model (for recent reviews, see Refs. pil| , pT| and references therein). In this paper, we 
investigate in a systematic way to what extent the high-energy part of the photon energy 
spectrum in B ^ Xg'-f decays can be controlled theoretically. The "Fermi motion" of 
the b quark inside the B meson, which determines the characteristic shape to the photon 
spectrum, can be consistently described by taking a convolution of the parton model 
prediction for the spectrum with a universal shape function F{k^), which determines the 



light-cone momentum distribution of the b quark in the B meson p6[-|]28[. We will for 
the first time present a discussion of Fermi motion effects in a full next-to-leading order 
analysis of B X^'y decays. In addition, our analysis contains several improvements 
over previous works concerning, in particular, the estimate of perturbative uncertainties, 
and the inclusion of QED corrections. In Section |^, we discuss in detail the structure of 
the B — > Xs'y branching ratio at next-to- leading order in QCD, correcting some errors 
in the formulae for real-gluon radiation contributions employed by previous authors. We 
introduce a decomposition of the branching ratio in terms of the values of the Wilson 
coefficients Ci{mw) at the weak scale, which is particularly convenient to discuss the 
sensitivity to New Physics beyond the Standard Model. This decomposition is also the 
starting point of a thorough discussion of the renormalization-scale dependence. We find 
that the perturbative uncertainty in the theoretical prediction for the branching ratio 
has been underestimated by previous authors by more than a factor of 2. We 

also suggest a new definition of the "total" branching ratio, which is insensitive to the 
unphysical soft-photon divergence in the b sg'-y subprocess. In Section ^, we show that 
Fermi motion effects, which result from the residual interaction of the b quark inside the B 
meson, give rise to the dominant theoretical uncertainty in the calculation of the partially 
integrated (over photon energy) branching ratio. We present a consistent treatment of 
these effects based on first principles of the heavy-quark expansion, emphasizing that for 
such partially integrated quantities the main element of uncertainty in the description of 
Fermi motion lies in the value of the 6-quark mass. Other features associated with the 
detailed functional form of the shape function play a minor role. We make a prediction 
for the B Xgj branching ratio with a restriction on the photon energy such that 
E}^^ > 2.2 GeV and find agreement with the CLEO measurement of the same quantity 
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within one standard deviation. We also extract a new value for the total branching 
ratio, which is significantly different from the published CLEO result reported in (|l|). In 
Section ^ we extend the discussion to the photon spectrum itself, investigating first the 
structure of the contributions from different operators in the effective Hamiltonian. We 
find that, to a high degree of accuracy, the shape of the photon spectrum is determined by 
QCD dynamics and is insensitive to New Physics beyond the Standard Model. We then 
perform a fit of our theoretical predictions to the CLEO data on the photon spectrum 
and extract again a value for the total branching ratio. We also discuss the possibility 
of determining, from future high-precision data on the photon spectrum, a value of 
the 6-quark mass with a well-defined short- distance interpretation. In Section ^ we 
consider the hadronic invariant mass spectrum and discuss the role of quark-hadron 
duality in the comparison of experimental data with our theoretical predictions. We 
derive a realistic, one-parameter description of the spectrum that is valid even in the 
low-mass region, where quark-hadron duality breaks down. Finally, in Section |^ we 
explore how New Physics beyond the Standard Model may affect the spectral shape 
and the total branching ratio in B ^ Xg'j decays. Section contains the conclusions. 
The paper also comprises three Appendices, where we discuss QED corrections to the 
B — »• Xgj branching ratio, the Doppler broadening of the photon spectrum in the decays 
of B mesons produced at the T(4s) resonance, and technical details of the calculation 
of the photon energy spectrum. 



2 B — > branching ratio 

The theoretical analysis of the B Xg'j branching ratio at next-to-leading order 
has been discussed previously by several authors. In this section we review the main 
ingredients of the calculation. In addition, we present several improvements of it, 
which concern the treatment of leading-logarithmic QED corrections, the analysis of the 
renormalization-scale dependence, and a discussion of the sensitivity to New Physics. 
We also correct some mistakes in the results for real-gluon emission presented in the 
literature. 

The starting point in the calculation of inclusive B decay rates is the low-energy 



effective Hamiltonian |^| 



-neflf = 



-^v,:v,,Y.c^{^bm^lb). (3) 

The operators relevant to our discussion are 

O2 = SL^f^CLCLl^bi , 



07 = ^ -SL^.uF'^'bn , 



O8 = f0^-L^..Gra«. (4) 
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To an excellent approximation, the contributions of other operators can be neglected. 
The renormalization scale fib in is conveniently chosen of order m^, so that all large 
logarithms reside in the Wilson coefficient functions. The complete theoretical prediction 
for the B —>■ Xsj decay rate at next-to-leading order has been presented for the first 
time by Chetyrkin et al. ||T^. It depends on a parameter 6 defined by the condition that 
the photon energy be above a threshold given hj > {1 — 5)E^^^, where E^^^ = mi,/ 2 
is the maximum photon energy attainable in the parton model. (Throughout this paper, 
we will neglect the mass of the strange quark whenever possible.) The prediction for 
the B X<j7 branching ratio is usually obtained by normalizing the result for the 
corresponding decay rate to that for the semileptonic decay rate, thereby eliminating a 
strong dependence on the 6-quark mass. We define 



T{B^Xsj)\ 



E-,>{l-S)E?^^^ 



Qa 



T{B ^ Xceiy) 



cb 



(5) 



where f{z) 



1 - 8^ + 8z^ 



0.542 - 2.23(v^ - 0.29) is a phase- 
space factor depending on the mass ratio z = {mc/mb)^, for which we shall take -y/i = 
0.29 ± 0.02. In the context of our analysis, the quark masses are defined as one-loop 
pole masses. The electro-magnetic coupling a = 1/137.036 is the fine-structure constant 
renormalized at = 0, as is appropriate for real-photon emission |Q. The quantity 
Ktsslo{^) = |C*7p + • • • contains the corrections to the leading-order result. In terms of 
the theoretically calculable ratio Rth{S), the B Xs'y branching ratio is given by 



B(5-.X,7)| 



E^>(l-S)E}} 



RttiS) X B{B ^X,eij) = O.lOSiVgL i?th(5) 



(6) 



where A^^sl = B(i? X^e u)/10.5% is a normalization factor to be determined from 
experiment. To good approximation A^sl = 1- The current experimental situation of 
measurements of the semileptonic branching ratio of B mesons and their theoretical 



interpretation are reviewed in Refs. ||3^, ^ . 

In the calculation of the quantity K-^i^oi^) '^^ shall consistently work to first order in 
the small parameters a^, 1/"^,^ and a/ag, the latter ratio being related to the leading- 
logarithmic QED corrections. The general structure of the result is 



i^NLo(5)= E %(5,/i,)ReCf)(/i,)Cf*(^,) +S{6) 



i,j=2,7,8 



+ SiS) 



2tt 



Re 



(0)*, 



a 



2 Re 



(0)*/ 



(7) 



where we have expanded the Wilson coefficients asQ 



Arc 



"^H/i.) + • • • • 



(8) 



^ These are the effective, scheme- independent Wilson coefficient functions introduced in Ref. |14| 
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The leading-order coefficients are given by 

("0), , 1 / _12 J>_\ 

i=i 

C^'\lx,)=V^^Ci'\mw) + j:hv'", (9) 

where 77 = as{mw)/ois{fib), and /ij, /ij and Oj are known numerical coefficients |T^, 0. 
In the Standard Model, the Wilson coefficients of the dipole operators at the scale mw 
are functions of the mass ratio Xt = {rnt{mw)/^wY given by |12 



C) \mw) = -jYz T^lnxt + ^ — '—^ 



4:{xt - 1)4 2A{xt - 1] 

The next-to-leading terms in @ must be kept only for the coefficient C-j{^h)- The expres- 
sion can be found in eq. (21) of Ref. [0. Our treatment of QED corrections 



differs from that of Czarnecki and Marciano in that we perform a renormalization- 
group improvement to resum the contributions to Cj{fib) of order aL (agL)^, with 
L = ln(mvK//^f)), to all orders in perturbation theory, whereas these authors include 
only the terms with n = 0. Numerically, the resummation decreases the effect of QED 
correction by almost a factor of 2. The technical details of our calculation are discussed 
in Appendix A. The result for (/if,) is 



^(em)/ X f o2 _9_ 4U _7_ 88 16\ (0). . 

( 32 _jL 32 _x 640 14 704 ^(0), 

-|- 77 23 -I ri 23 -I 7723 -0 23 Co mw 

V 575 ' 1449 ' 1449 ' 1725 ' / « v 



190 „35 359 _i7 4276 _i2 350531 _^ 

■ 7y 23 — 7^ 23 -] 7^ 23 -| 23 



8073 ' 3105 ' 121095 ' 1009125 

2 _^ 5956 A 38380 i4 748 le 

-| fi 2a — 7723 -I- 7123 — 7723 

4347 ' 15525 ' 169533 ' 8625 ' 



The result for the leading QED correction to the semileptonic decay rate is 



(11) 



4r(/^.) = ^ [V-' - 1 = In ^ . (12) 

23 ^ ^ TT yUf, 
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We note that unlike the factor a in (|^), which results from the calculation of a matrix 
element for a process with real-photon emission, the QED corrections to the Wilson 
coefficients arise from the evolution of local operators, and hence the coupling a in (^) 
should in principle be taken as a running coupling rather than the fine-structure 
constant renormalized at = 0. However, including the running of the QED coupling in 
the operator evolution would only induce corrections of order (aL)^(asL)", which from 
a numerical point of view can be safely neglected. 

For the purpose of illustration, we note that with /if, = 4.8 GeV the values of the 
various coefficients in the Standard Model are: C2\mh) ~ 1.11, cf^imb) ~ —0.31, 
Cf\mb) ^ -0.15, as well as Ci^\mb) « 0.48 and Ci'"^\mb) ^ 0.03. The QED correction 
proportional to Cj^"^^ in is about a factor 0.13 smaller than the next-to-leading order 
QCD correction proportional to Cj^\ 

The coefficient functions kij{6,fib) in are given by 



kniS, fib) 




rj + In 



+ 



1 



6A2 



f77{S)+S{6) 



27T 



k27iS, /ife) 

kij{6, Hb) 



8(5) 



SiS) 



27r 

27T 









A2 ■ 


f^bj 




f^bj 


(ysifJ'b) 


71 



H l27[d) 



7T 



C^sifJ'b) 



71 



where 



S{5) = exp 



f7s{S) 

{z,j} = {2,2}, {8,8}, {2,8}, 
2aM ^1,2^ + 1 1,5^ 



is a Sudakov factor, 777 = 
dimension matrix, and 

Re(r2) 



727 



10 

y 



416 
81 



37T 

and 787 



(13) 



(14) 



Y are entries of the anomalous 



Re(r8) 



44 



87r2 
27 



-4.092 + 12. 78(v^- 0.29) 



(15) 



are numerical coefficients resulting from the calculation of the matrix elements of the 
local operators Oi in the effective Hamiltonian at next-to-leading order [0. Finally, 
R,{z) ~ 3.382 — 4.14(yi — 0.29) is the next-to-leading correction to the semileptonic 
decay rate To 0{as) the explicit fib dependence of the coefficients kij{S,fib) cancels 
against that of the Wilson coefficient functions. Following Ref. we allow for different 
renormalization scales in radiative and semileptonic B decays (i.e. fib 7^ fib)- 
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The functions fij{S) in (|T3|) account for the effects of real-gluon radiation and are 



defined such that fij{0) = 0. They can be obtained from results given in Refs. [|T7[ [18 
by performing some phase-space integrations. We find 



MS) 



/78(5) 
/22(5) 



105 + - — + 6(5 - A) \n5 



^ Ul2{1 -6) -^ + 8 hi(l -6)- 6{2 + 6)\n6 



+ 76 + 35^ 



L2{l-6) 



25'^ 



25 + 5^ + A\n{l - 5) 



In 



51n5H ^ 

6 4 4 12 



16 
27 



dx (1 — a;)(l — x^) 



X \zj 2 



/27(5) = -3/28(5) = -y /da; (1 - x^) Re 



(16) 



where X5 = max(x, 1 — 5), and 



G{t) 



-2arctan%/t/(4 - t) 



2 In 



{Vi+Vt^)/2 



t < 4, 



t > 4. 



(17) 



Our expressions for fjsi^) and fssi^) disagree with the corresponding ones in Ref. [|T^ , 
which have later been used by several authors. (The corrected expressions are also given 
in an Erratum to Ref. [0.) We shall comment on the numerical effect of this correction 
below. The function /gs is sensitive to collinear singularities regulated by the mass of the 
strange quark. The collinear logarithms can be resummed to all orders of perturbation 
theory, leading to a coUinear-safe result pSf. Unless 5 is chosen very close to 1, the net 



effect of the resummation is a moderate increase of the result. Since the contribution 
proportional to fss is very small, however, it is sufficient for all practical purposes to 
work with the leading-order expression given above. We take a rather large value for 
the quark-mass ratio, mb/rris = 50, in order to mimic the effect of the resummation of 
collinear logarithms. 

Bound-state corrections enter the theoretical expressions for the coefficients kij at 
order l/mq and are proportional to the hadronic parameter A2 = ^i^ri^, — m\) ^ 
0.12GeV^ 1^91 . The corrections proportional to 1/ml entering the expression for kn- 
characterize a spin-dependent interaction inside the B meson A peculiar feature of 

inclusive radiative decays is the appearance of the correction proportional to 1/m^ in /C27, 
which represents a long-distance contribution arising from (cc) intermediate states 0, |^. 
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Strictly speaking, these effects are non-local in nature; however, to a good approximation 
they can be represented by a local correction [1^. The correct sign of this term 



has only recently been found in Ref. [11 



Finally, a comment is in order about our treatment of the Sudakov factor, which is 
slightly different from that of previous authors. In Refs. |19[-||21|, the Sudakov factor was 
only included for the leading term in fcyy. If this is done, the decay rate becomes negative 
for small values of 6, which is an unphysical result. All terms not vanishing in the limit 
6 = correspond to a two-body decay b ^ S'j and must be suppressed by a Sudakov 
factor. Once this is done, the quantity ii^NLol*^) vanishes in the limit 5 — as it should. A 
full next-to-leading order resummation of Sudakov logarithms that goes beyond the naive 
exponentiation of the one-loop result shown in (|I^ is possible but rather complicated. 
In Refs. EDI, ET|, such a resummation has been performed for high-order moments of the 



photon energy spectrum; however, the results are such that a numerical evaluation would 
require integration over the running coupling constant as{k±) in the region Aqcd < fc± < 
rUb- In Ref. the resummation has been extended to the partially integrated photon 
spectrum itself (rather than its moments), and a factorization of short- and long-distance 
contributions has been performed such that all contributions from scales k_i_ in the range 
(^QCD^fc)^^^ < kj_ < rrih are treated perturbatively, whereas contributions from scales 
in the range Aqcd < k± < {AqcDf^bY^'^ are absorbed into the definition of the shape 
function. This guarantees that the resummed formulae can be reliably evaluated in 
perturbation theory. If this is done, it turns out that the resummation is a very small 
effect, which can be neglected for all practical purposes [Q. 

In order to explore the sensitivity of the theoretical prediction for the B Xg'y 
branching ratio to possible New Physics contributions, it is instructive to make explicit 
the dependence of the result on the values of the Wilson coefficients of the dipole oper- 
ators O7 and Og at the scale m^y. To this end, we introduce the ratios 

They are normalized to the Standard Model contributions, which at next-to-leading 
order take the values Cf^(mH/) ~ —0.22 and C^^{mw) ~ —0.12 [0, |I6|. In many 



extensions of the Standard Model there are contributions to Cj{mw) and C^i^mw) from 
new flavour physics at a high scale, and consequently the parameters 1^7 and may 
take (even complex) values different from 1. Similarly, New Physics may induce dipole 
operators with opposite chirality to that of the Standard Model, i.e. operators with 
right-handed light-quark fields. If we denote by and the Wilson coefficients of 
these new operators, the expression (^ can be modified to include their contributions 
by simply replacing CjC* ^i^j + C^C^* everywhere, taking however into account 
that C2 = 0. We thus define two additional parameters 

which vanish in the Standard Model. Since the dipole operators only contribute to rare 
fiavour- changing neutral current processes, there are at present rather weak constraints 
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on the values of these parameters (see Section On the other hand, we assume that 
the coefficient C2 of the current-current operator O2 takes its Standard Model value, 
and that there is no similar operator containing right-handed quark fields. Since the 
operator O2 mediates Cabibbo-allowed decays of B mesons, any significant New Physics 
contribution to C2 would already have been detected experimentally. 

With these definitions, the B —>■ Xs'j branching ratio can be decomposed as 

B{B Xs'j)\e^>(i-5)e;:^^^ 
= B22{5) + B,,{5) + + B,,{5) (l^sp + 

+ B2-j{5) Re(e7) + B2s{5) Re{i^) + B^^) \^<i7il) + ^6(^7%''*)] • (20) 

The components Bij{6) are formally independent of the renormalization scale fib- Their 
residual scale dependence results only from the truncation of perturbation theory at next- 
to-leading order. In Table |^, the values of these quantities are given for different choices 
of the renormalization scale and the cutoff on the photon energy. The input parameters 
entering the calculation will be discussed below. Typically, the components Bij vary by 
amounts of order 10-20% as /if, varies between mb/2 and 2m;,. The good stability is a 
result of the explicit cancelation of the fib dependence between the Wilson coefficients 
and matrix elements achieved by a full next-to-leading order calculation. The Standard 
Model branching ratio is obtained by adding the various contributions setting = C,s = ^ 
and = ~ 0; shown in the last column. The most important contributions are 
the 2-2 and 2-7 terms, followed by the 7-7 term. Note that with a realistic choice of the 
cutoff parameter 6 the coefficient B^siS) of the term proportional to l^gp + is very 
small. Therefore, B —>■ Xg'j decays have a low sensitivity to enhanced chromo-magnetic 
dipole transitions. For the remainder of this section we focus on the Standard Model and 
evaluate the various theoretical uncertainties in the prediction for the branching ratio. 
The impact of New Physics will be discussed in Section ^. 

In Table |1|, the choice 6 = 0.9 corresponds to the unrealistic case of an almost 
fully inclusive measurement, whereas 6 = 0.3 and 0.15 correspond to restrictions to the 
high-energy part of the photon spectrum, which in practice is required for experimental 
reasons. The theoretical prediction for the branching ratio diverges in the limit 6 1 
because of a logarithmic singularity in the term proportional to fssi^), which reflects 
the soft-photon divergence of the b — >■ sgj subprocess. Previous authors |T^-|21[] have 



chosen to define the "total" B — > Xg'y branching ratio by taking 6 = 0.99. In our opinion 
this is not the best definition possible, because for values of 6 so close to 1 the theoretical 
result becomes very sensitive to the unphysical soft-photon divergence. This is evident 
from Figure m, which shows the integrated branching ratio as a function of 6. We believe 
a more reasonable definition of the "total" branching ratio is to use an extrapolation to 
6 = 1 starting from the region 6 ~ 0.5-0.8, where the theoretical result exhibits a weak, 
almost linear dependence on the cutoffi The simple geometric construction indicated by 
the dashed lines in the figure shows that the extrapolated value so defined agrees, to a 
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Table 1: Values of the coefficients Bij{6) in units of 10 ^, for different choices of fi^ 



fib 5 


B22 B-j-j Bs8 B27 B28 -B78 


j:b., 


mb/2 0.90 
0.30 
0.15 


1.321 0.335 0.015 1.265 0.179 0.074 
1.167 0.322 0.005 1.196 0.136 0.070 
1.080 0.309 0.004 1.143 0.126 0.067 


3.188 
2.896 
2.728 


rub 0.90 
0.30 
0.15 


1.258 0.382 0.015 1.395 0.161 0.083 
1.239 0.361 0.005 1.387 0.137 0.080 
1.200 0.347 0.004 1.354 0.132 0.077 


3.293 
3.210 
3.114 


2mb 0.90 
0.30 
0.15 


1.023 0.428 0.015 1.517 0.132 0.092 
1.041 0.402 0.004 1.552 0.118 0.091 
1.021 0.386 0.004 1.535 0.115 0.088 


3.206 
3.209 
3.150 



3.5 



I-:) 




s 

Figure 1: Dependence of the B Xg'j branching ratio on the cutoff parameter 6. The 
dashed lines indicate the extrapolation to the "total" branching ratio. 



good approximation, with the result obtained by taking 5 = 0.9. Hence, from now on 
we define the "total" branching ratio to be that corresponding to this particular value 
of the cutoff. 

The dependence of the theoretical results on the choice of the renormalization scale 
is conventionally taken as an estimate of higher-order corrections. Following common 
practice, we vary the renormalization scales fib and fib independently in the range between 
mb/2 and 2mb; their central values are taken to be nib- The dependence on the scale 
fib entering the formula for the semileptonic decay rate is straightforward to analyse. 
The result is shown in Table |^. The analysis of the scale dependence of the radiative 
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decay rate is more subtle. Previous authors have estimated the jj^b dependence of the 
total branching ratio in the Standard Model and found a striking improvement over the 
leading order result. A dedicated discussion of this issue has been presented by Buras 



et al. pOl' results obtained in this work were later confirmed in Refs. |21|, p^ . 

One finds a variation of the total branching ratio by l3;2%, as compared with at 
leading order. We agree with those results, but we believe they cannot be taken as a 
realistic estimate of the size of unknown higher-order corrections. The excellent stability 
observed at next-to-leading order is largely due to an accidental cancelation between 
different contributions to the decay rate. This point of view is supported by the fact 
that in some extensions of the Standard Model, such as two-Higgs-doublet models, a 
much stronger scale dependence is observed in some regions of parameter space . 
On the left-hand side in Figure 0, we show the branching ratio as a function of fib/rrih, 
both at leading and at next-to-leading order. The leading-order result is obtained by 
replacing the quantity -f^NLol*^) with \c!f\jjb)\'^ . The three curves in the upper plot 
refer to different choices of 6 (there is no cutoff dependence of the leading-order result). 
Note the different scales in the two plots. The improvement in going from the leading 
to the next-to-leading order is spectacular and reduces the apparent scale dependence 
by more than a factor of 10. However, a more careful look at the upper plot reveals 
two surprises: first, the scale dependence increases rapidly as /if, is taken below 0.7m;,, 
although perturbation theory should work well for lower scales than that; secondly, 
for fib > 2mb the prediction for the partially integrated branching ratio with 5 = 0.3 
exceeds the prediction for the total branching ratio (obtained with 6 = 0.9), which is an 
unphysical result. Both observations indicate that higher-order corrections may be more 
important than what is suggested by the apparent weak scale dependence of the curves 
in the plateau region. The scale dependence of the three most important contributions 
to the branching ratio (those from B22, -B27 and -B77) is illustrated in the right-hand 
plots in Figure There is again a significant improvement in going from the leading 
to the next-to-leading order. However, the residual scale dependence of the quantities 
Bij at next-to-leading order is much larger than that of their sum, which determines 
the total branching ratio in the Standard Model. Note, in particular, the almost perfect 
cancelation of the scale dependence between the 2-2 and the 2-7 terms, which is accidental 
since the magnitude of the 2-7 term depends on the top-quark mass through the value 
of Ci'^ {mw) in (|10D, whereas the 2-2 term is independent of m^. In such a situation, the 
apparent weak scale dependence of the sum of all contributions is not a good measure 
of higher-order corrections. Indeed, higher-order corrections must stabilize the different 
curves in the right-hand upper plot individually, not only their sum. The variation of 
the individual components Bij as a function of fib thus provides a more conservative 
estimate of the truncation error than does the variation of the total branching ratio. For 
each component, we estimate the truncation error by taking one half of the maximum 
variation obtained by varying fib between mb/2 and 2mb. The truncation error of the 
sum is then obtained by adding the individual errors in quadrature. As shown in Table ^ 
we find a total truncation error of about ±7%, which is more than a factor of 2 larger 
than the estimates obtained by previous authors |T3-P^. In our opinion an even larger 
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Figure 2: Scale dependence of the B Xs'j branching ratio (left) and of its three 
most important components (right) in the Standard Model. For comparison, the lower 
plots show the results obtained at leading order. 



truncation error could be justified given that the choice of the range of variation of /ib is 
ad hoc, and that the scale dependence of the various curves in Figure ^ is not symmetric 
around the point fib = "^6- 

Let us finally discuss the sensitivity of the theoretical prediction for the B Xs'j 
branching ratio to the various input parameters entering the calculation. For the quark 



Table 2: Different sources of theoretical uncertainties (in 



6 



0.90 
0.30 
0.15 



fib fnc/mb nib rut as{mz) CKM EW cor. 



±6.3 
±6.6 
±7.1 



+2.2 


+5.9 


-1.0 


+1.6 


+2.8 


-1.5 


-5.0 


+1.1 


-1.7 


-2.7 


+2.2 


+5.5 


-1.0 


+1.7 


+2.6 


-1.5 


-4.6 


+1.0 


-1.7 


-2.6 


+2.2 


+5.3 


-0.9 


+1.6 


+2.5 


-1.5 


-4.5 


+1.0 


-1.7 


-2.4 



±2.1 ±2.0 
±2.1 ±2.0 
±2.1 ±2.0 



total 



-10.0 

- 9.3 
f9.9 
-9.3 
-10.1 

- 9.5 
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pole masses, we take -J~z = nic/mb = 0.29±0.02, nib = (4.80±0.15) GeV, and mt = (175± 
6) GeV corresponding to the running mass mt{mw) = (178±6) GeV. We use the two-loop 
expression for the running coupling Q;s(/i) with the initial value as{mz) = 0.118 ±0.003. 
Finally, for the ratio of the CKM parameters we take the value | V^^Vtbl/l Vcfe| = 0.976 ± 
0.010 obtained from a global analysis of the unitarity triangle |^^. We also include 
an uncertainty of ±2% to account for next-to-leading electroweak radiative corrections 
||33| . The theoretical uncertainties arising from the variation of each input parameter are 
collected in Table ^ Adding the different errors in quadrature we get total uncertainties 
of about ±10% in all three cases. For the total branching ratio in the Standard Model 
we obtain B{B Xg'y) = (3.29 ±0.33) x 10~^iVsL. Contrary to common folklore, about 
40% of the branching ratio reflects the presence of the current-current operator O2 in 
the effective Hamiltonian and is not related to penguin diagrams with a top-quark loop. 
For comparison with previous authors, we note that with the choice 6 = 0.99 we would 
obtain B{B X,7) = (3.37 ± 0.34) x IO-^A^sl. 



3 Implementation of Fermi motion 

Whereas the explicit power corrections included in the functions ku and in (0) 
are very small, there is an important nonperturbative effect that has not been included 
so far in next-to-leading order analyses of the B Xs'j branching ratio: the residual 
motion of the b quark inside the B meson caused by its soft interactions with the light 
constituents leads to a modification of the photon energy spectrum, which is an impor- 
tant effect if a realistic cutoff is imposed |0. This so-called "Fermi motion" is included 



in the heavy-quark expansion by resumming an infinite set of leading-twist corrections 
into a shape function F(/c+), which governs the light-cone momentum distribution of 
the heavy quark inside the B meson |26[-|^. This function shares many similarities 



with the parton distributions in deeply inelastic scattering. The physical decay distri- 
butions are obtained from a convolution of parton model spectra with this function. In 
the process, phase-space boundaries defined by parton kinematics are transformed into 
the proper physical boundaries defined by hadron kinematics. The shape function is a 
universal, i.e. process-independent characteristic of the B meson governing the inclusive 
decay spectra in processes with massless partons in the final state, such as B ^ Xs'j and 
B — > Xu i V- It is important to note that this function does not describe in an accurate 
way the distributions in decays into massive partons such as i? — > Xcd-v p7| , [28[] . Unfor- 
tunately, therefore, the shape function cannot be determined using the lepton spectrum 
in semileptonic decays of B mesons, for which high-precision data exist. On the other 
hand, there is some useful theoretical information on the moments of the shape function, 
which are related to the forward matrix elements of local operators : 



An = I dk+ kl F{k+) = {B\ b {iD+Yb \B) . (21) 



2m 



B 



The first three moments satisfy Aq = 1, Ai = and A2 = where //^ = — Ai is 
related to the kinetic energy of the b quark inside the B meson [^. The condition 
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Ai = 0, which is a consequence of the equations of motion, ensures that the quark mass 
rrih entering our theoretical expressions is the pole mass (defined to the appropriate order 
in perturbation theory, i.e. to one- loop order for our purposes). 

Let PpiVp) be the photon energy spectrum in the parton model, where Up = 2E^/mb 
with < i/p < 1. Our goal is to include the effects of Fermi motion and calculate the 
physical spectrum P{y) as a function of the variable y = 2E^/mB. To leading-twist 
approximation, the result is given by the convolution |^ 



P{y)dy=[dk+F{k+)\PMdyp\ _ ,,,, (22) 
where yp{k^) is obtained by replacing mb in the definition of yp with the "effective mass" 



ml = mb + i.e. |/p(fc+) = 2E^/ml = yniB/ml. Because the support of the shape 

function is restricted to the range —irtb < /c+ < — m;,, it follows that < y < 1. 
In other words, after the inclusion of Fermi motion the spectrum extends to the true 
kinematic endpoint at E^^^ = itib/^. Let us denote by i?p(5p) the integrated branching 
ratio calculated in the parton model, which is given by an integral over the spectrum 
Pp{yp) with a cutoff 5p defined by the condition that E^^ > |(1 — 5p)mb. From (|2^), 
it then follows that the corresponding physical quantity B{6) with 6 defined such that 
> |(1 — S)'mB is given by 

BW= / ,k,Fik,)B,[,-^^lsLzJ.y (23) 

ms(l-<5)-»ni, 

This relation is such that -B(l) = -Bp(l), implying that the total branching ratio is 
not affected by Fermi motion; indeed, the I/itlq corrections in (|T3[) are the only power 
corrections to the total branching ratio. The effects of Fermi motion are, however, 
important for realistic values of the cutoff. We will now evaluate relation ( P5| ) for the 
various components Bij{6) introduced in the previous section. 

Several ansatze for the shape function have been suggested in the literature [^ - [^ . 
For our purposes, given the poor present knowledge about higher moments An with 
n > 3, it is sufficient to adopt the simple form 

F{k+) = N{1- a;)'^e(^+'^)^ ; a; = ^ < 1 , (24) 

where A = tub — mi,. This ansatz is such that Ai = by construction,^ whereas the 
condition Aq = 1 fixes the normalization A^. The parameter a can be related to the 
second moment, yielding A2 = |/U^ = A^/(l + a). Thus, the 6-quark mass (or A) and 
the quantity f^l (or a) are the two parameters of our function. Below, we will take 
mb = 4.8 GeV and /x^ = 0.3 GeV^ as reference values, in which case a ~ 1.29. Because 
the main effect of Fermi motion is to fill the gap between the parton model endpoint 
of the photon spectrum and the physical endpoint, it turns out that the results are 

•^For simplicity, we neglect exponentially small terms in nib/ A. 
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Table 3: Values of the coefRcients Bij{6) in units of 10 corrected for Fermi motion 



mb [GeV] 6 Ef"" [GeV] 


B22 Bj7 Bss B27 B28 Bjs 


EBij 


4.65 0.90 0.26 
0.30 1.85 
0.15 2.24 


1.289 0.375 0.014 1.401 0.162 0.083 
1.201 0.333 0.004 1.322 0.132 0.075 
0.802 0.220 0.002 0.889 0.087 0.050 


3.324 
3.068 
2.050 


4.80 0.90 0.26 
0.30 1.85 
0.15 2.24 


1.258 0.382 0.014 1.395 0.160 0.083 
1.213 0.352 0.004 1.362 0.134 0.078 
0.945 0.272 0.003 1.071 0.104 0.060 


3.291 
3.144 
2.456 


4.95 0.90 0.26 
0.30 1.85 
0.15 2.24 


1.227 0.388 0.014 1.388 0.157 0.083 
1.200 0.365 0.004 1.375 0.133 0.080 
1.072 0.323 0.003 1.239 0.118 0.070 


3.259 
3.156 
2.825 



very sensitive to the choice of the 6-quark mass. Table |^ shows the coefficients Bij{S) 
corrected for Fermi motion using the above ansatz with a fixed value a ~ 1.29 but 
different values of m^. A representative range of parameters is covered by considering 
the three following cases: mt = 4.65 GeV (yielding A ^ 0.63 GeV and /i^ ^ 0.52 GeV^), 
nib = 4.8 GeV (yielding A ^ 0.48 GeV and fil ^ 0.3 GeV^), and rrib = 4.95 GeV (yielding 
A ^ 0.33 GeV and /i^ ~ 0.14 GeV^). We also show the values of the photon energy 
cutoff, -E™™ = |(1 — 6)mB, which are now independent of the 6-quark mass. 

For a graphical illustration of the sensitivity of our results to the parameters of the 
shape function, we show in the upper plots in Figure ^ the predictions for the Standard 
Model branching ratio as a function of the energy cutoff i?™". In the first plot, we 
use the same sets of parameters as in Table |, i.e. rrib = 4.65 GeV (long-dashed curve), 
4.8 GeV (solid curve), and 4.95 GeV (short-dashed curve), with /i^ adjusted such that the 
ratio /i^/A^ remains constant. The gray line shows the result obtained using the same 
parameters as for the solid line, but with a Gaussian ansatz F{k^) = N {1 — x)'^e~^^^~^'' 
for the shape function. For comparison, we also show the data point B(S — Xs'~f) = 
(2.04 ± 0.47) X 10"^ obtained by the CLEO Collaboration with a cutoff at 2.2 GeV [|3 



As shown in Appendix B, the fact that in the CLEO analysis the cutoff is imposed 
on the photon energy in the laboratory frame rather than in the rest frame of the B 
meson is not very important for the partially integrated branching ratio and has been 
neglected here. In the second plot, we keep = 4.8 GeV fixed and compare the parton 
model result (gray dotted curve) with the results corrected for Fermi motion, using 
nl = 0.15 GeV^ (short-dashed curve), 0.30 GeV^ (sohd curve), and 0.45 GeV^ (long- 
dashed curve). This figure illustrates how Fermi motion fills the gap between the parton 
model endpoint at mb/2 and the physical endpoint at 1713/2. To be precise, the physical 
endpoint is actually located at [m^ — (m^ + m.j,Y]/2mB ~ 2.60 GeV, i.e. slightly below 
ms/2 ^ 2.64 GeV. Close to the endpoint, our theoretical prediction provides an average 
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Figure 3: Theoretical predictions for the integrated B — > Xg^f branching ratio (upper 
plots) and the corresponding photon spectra (lower plots) for various choices of the 
shape-function parameters (mh, ^\) and functional form, as explained in the text. The 
calculation of the photon spectra will he discussed in Section ^. 



description of the true spectrum in the sense of quark-hadron duahty (see Section |^). 
Comparing the two upper plots in Figure we observe that the uncertainty due to 
the value of the 6-quark mass is the dominant one. Variations of the parameter /i^ 
have a much smaller effect on the partially integrated branching ratio, and also the 
sensitivity to the functional form adopted for the shape function turns out to be small. 
This behaviour is a consequence of global quark-hadron duality, which ensures that even 
partially integrated quantities are rather insensitive to bound-state effects. The strong 
remaining dependence on the fe-quark mass is simply due to the transformation by Fermi 
motion of phase-space boundaries from parton to hadron kinematics. We believe that the 
spread of results obtained by varying rrif, between 4.65 and 4.95 GeV (with fi^ adjusted 
as described above) is a fair representation of the amount of model dependence resulting 
from the inclusion of Fermi motion. With a cutoff -E™™ = 2.2 GeV as used in the CLEO 
analysis, and correcting for the small effect of the boost from the B rest frame to the 
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laboratory frame (see Appendix B), we obtain 

B(fi ^ X,7)U.a.>2.2GeV = (2-57 ± 0.26t°:i) x IQ-^iVgL . (25) 

The first error accounts for the dependence on the various input parameters, while 
the second one reflects the uncertainty due to the modeling of Fermi motion. For 
a cutoff as high as that employed in the CLEO analysis, this uncertainty is in fact 
the dominant theoretical error. In the future, an effort should therefore be made to 
lower the cutoff on the photon energy to a value of 2 GeV or less. Comparing our 
result with the CLEO measurement of (2.04 ± 0.47) x 10~^ we obtain the ra- 
tio R = Bexp/Bth = 0.79 ± 0.18(exp)lo;j2(th) (assuming A^^sl = 1), which deviates 
from unity by less than one standard deviation. This must be confronted with the 
comparison of the total branching ratios using the value reported in (|^), which gives 
R = Boxp/Bth = 0.71 ± 0.20(exp) ± 0.07(th). The extrapolation to low photon ener- 
gies performed in the CLEO analysis [|l| has artificially increased the deviation from the 
theoretical prediction by a significant amount, i.e. the model dependence inherent in 
this extrapolation has been underestimated. This is also reflected in the fact that the 
correction factor relating the total branching ratio to the branching ratio obtained with 
the cut E\^^ > 2.2 GeV is K2.2 = 0.78 

-ail! which differs significantly from the factor 
0.87 ± 0.06 employed by CLEO. Using our correction factor, the extrapolation of the 
CLEO measurement to the total branching ratio yields 

B{B ^ X,7)cLEO = (2.62 ± 0.60exp ^[JJ^th) x 10"^ (26) 

instead of the value quoted in ([^). We stress that the change in the central value and 
the increase of the theoretical errorQ with respect to the result reported by CLEO are 
entirely due to the improved treatment of bound-state effects presented in this paper. 
Whereas the CLEO analysis relies on a quark model, we perform an analysis that is 
entirely based on QCD and the operator product expansion. Our treatment is thus not 
only more conservative but also more consistent from a theoretical point of view. 

The procedure of extrapolating a measurement of the B —>■ Xs'~f branching ratio in 
the region of high photon energies to the total branching ratio not only introduces large 
systematic errors, but also entails the disadvantage that one has to rely on the Standard 
Model to describe the photon spectrum in the low-energy region. In our opinion, it 
would therefore be desirable if in the future the comparison of theory with experiment 
were done for the partially integrated branching ratio, which is the quantity actually 
measured, rather than for the total branching ratio. 

4 Photon spectrum and determination of 

The large theoretical uncertainties in the calculation of Fermi motion effects on the 
partially integrated branching ratio for B Xg'y decays can be reduced in two ways. 

"^The second error quoted in (0) is dominated by experimental systematics. The theoretical error 
was assumed to be ±0.16 x 10^*. 
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The first possibility is to lower the cutoff i?™™ on the photon energy. As is apparent 
from Figure if a value as low as i?™™ ^ 2 GeV could be achieved, the theoretical 
predictions would become insensitive to the parameters of the shape function. To what 
extent this will be possible in future experiments depends on their capability to reject the 
background of photons from other decays. The Cabibbo-favoured B decays into charmed 
particles, in particular, can yield photons of energy up to about 2.3 GeV. The second 
possibility is that future high-precision measurements of the photon spectrum will make it 
possible to adjust the parameters of the shape function from a fit to the data. We repeat 
that these parameters cannot be determined from a study of the lepton spectrum in 
B — > Xc i V decays. On the other hand, a determination of the shape-function parameters 
from B — > XsPi decays would enable us to predict the lepton spectrum in 5 iv m. 

a model- independent way [^. This may help to reduce the theoretical uncertainty in 
the current value of A detailed analysis of the photon spectrum will therefore be 

an important aspect in future analyses of inclusive radiative B decays. 

Given the expression for the integrated B X^.'^ branching ratio derived in the 
previous sections, the photon spectrum can be obtained from differentiation with respect 
to 5, i.e. 



1 



B(5 ^X^ev) 



dB(5 X,7) _ 6a 

dy 'Kf{z) 



ch 



i^^,Lo(l-Z/)> (27) 



where y 



E^/E^^^. In analogy with (^, we write 

NLo(l -?/) 



a 



i,j=2,7,8 



271 



Re 



(0)* 



a 



2 Re 



(0)*/ 



where 
A(y) 



( ln(l - y) + 7 ) exp 



2aa(/ifc 
37" 



y) 



(28) 



(29) 



is the derivative of the Sudakov factor. The Sudakov factor can be regarded as a smeared 
step function, and hence A(?/) can be viewed as a smeared (5-function. In other words, we 
must consider A(?/) = 0(1) rather than 0{as)-, in spite of the prefactor in (^9]). The 
coefficient functions Pij{y , jJ'b) in (^ 



have the same form as the coefficients kij{S, /j^b) 



m 

f 



Tdj ) but with the replacements S{6) — >• A(?/) and fij{S) — > Sij{y), where Sij{y) 



jj^^ — y). The explicit expressions for these functions are presented in Appendix C. 
Given the photon spectrum Ppiyp) in the parton model, where yp = 2Ey/mb, the 
next step is to implement Fermi motion. According to (p2l), this is achieved by taking 
the convolution 



niB-mi, 



P{y) 



rriB 



niBy-mb 



nib + ^4 



msy 
nib + 



(30) 
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where y = 2E^/mB- It is convenient to decompose the final result for the spectrum in a 
form analogous to (^) by writing 

1 dB(B^A-.oO_oio5^Ap(2£,/,„,) 



-R\2 



NsL dE^ niB 

= P,,{E,) + P,,{E,) [U' + + P,s{E,) + \ii 

+ P2^{E,) Re(e7) + ^28(^^7) Re(e8) + ^78(^^7) [^^(^7^) + Re(e7%''*)] . (31) 

The results for the various components of the spectrum are shown in Figure ^, where 
we take central values of all input parameters. The contributions are ordered according 
to their magnitude. In the last plot, we show all components together on a logarithmic 
scale. Note that, with the exception of the tiny 8-8 contribution, the different components 
have a very similar spectral shape. This observation implies that the shape of the photon 
spectrum is not sensitive to physics beyond the Standard Model. With a realistic cutoff 
on the photon energy, even large deviations of the parameters ^7^'' and ^g^^ from their 
standard values would not have a detectable effect on the shape of the photon spectrum. 
Although this may be disappointing from the point of view of searching for New Physics 
in 5 — > Xs'^ decays, it entails the advantage that a precise measurement of the spectrum 
can be used to determine the parameters of the shape function without relying on the 
Standard Model. For the remainder of this section we concentrate on the Standard 
Model, for which the photon spectrum is given by the sum of the individual contributions 
shown in Figure ^ The results obtained for various choices of the parameters of the shape 
function are shown in the lower plots in Figure ^. The photon spectra are more sensitive 
to the functional form of the shape function than are the predictions for the integrated 
branching ratio in the upper plots. Therefore, a fit to future high-precision data on the 
spectrum should use a more fiexible ansatz for the shape function than the one given 
in (p^. On the other hand, we will see below that even a small element of smearing 
provided, e.g., by the finite detector resolution or the Lorentz boost of photons from the 
B rest frame to the laboratory frame, is sufficient to reduce this sensitivity significantly. 

Let us now make a comparison of our predictions for the photon energy spectrum 
with the data obtained by the CLEO Collaboration, which are presented in Table ^ 
To this end, we must take into account that B mesons produced in decays of the T(4s) 
resonance have a small momentum in the laboratory system, so that the photon spectrum 
is Doppler shifted. The boost that connects the B rest frame with the laboratory frame 
is characterized by 



IpbI 

En 



\ 



^ = 0.064 ±0.007, (32) 



2 

"^T{4s) 



and the maximum energy in the laboratory frame is {E}^^)^^^,^ ^ {1 + (3)E^^^ ^ jTjrnax _|_ 
170 MeV. The formalism for incorporating this effect is discussed in Appendix B. The re- 
sults are shown in the left-hand plot in Figure |^, where we compare the corrected photon 
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Figure 4: Different components of tlie photon spectrum in B Xg^ decays 



spectra with the CLEO data, using our standard parameters for the shape function. No 
fit to the data has been performed; in all cases, the theoretical spectra are normalized to 
the central Standard Model value of 3.29 x 10~^ for the total branching ratio (assuming 
A'^sL = 1). Note that after the smearing implied by the Doppler shift of the spectra, 
the gray and the black solid lines, which as before correspond to different functional 
forms adopted for the shape function, are very close together. This reflects the reduced 
sensitivity to the fine details of the modeling of Fermi motion, which is achieved by any 
kind of smearing of the photon spectrum. 

To perform a fit to the data, we rebin our theoretical results in the same energy 
intervals as used by CLEO and, for each set of parameters for the shape function, adjust 
the overall normalization (i.e. the total branching ratio) to give the best fit to the data. 
The results are reported in Table ^, and the best fits displayed in the right-hand plot 
in Figure ^ All fits have an excellent x^/'^dof ^ 1? indicating that with the present 
accuracy of the data it is not possible to determine the parameters of the shape function 
in a meaningful way. Unless a very large value of m?, is chosen, the result for the total 
branching ratio comes out higher than the value (|l]) reported by CLEO, and the upper 



20 



Table 4: CLEO results for the photon spectrum in B — > X^'^ decays p^/ . The value 
of the branching ratio reported in ^ was obtained from the middle two bins. 





[GeV] 


dB/d^^^^ [10-^GeV"^] 


1.95- 


-2.20 


2.13 ±2.38 


2.20- 


-2.45 


4.50 ± 1.52 


2.45- 


-2.70 


3.54 ±0.98 


2.70- 


-2.95 


0.11 ±0.54 




2.2 2.4 2.6 2.8 
E}f' [GeV] 



> 
O 
I 

o 



mt [GeV] : 

4.65 

4.80 

4.95 



2.2 2.4 2.6 
^lab 



[GeV] 



Figure 5: Theoretical predictions for the photon energy spectrum in the laboratory 
frame for different parameters of the shape function. The gray line in the left-hand 
plot shows the result obtained using a Gaussian form of the shape function with central 
values of mi, and /U^. The data points show the CLEO results. In the left-hand plot, 
no ht to the data is performed, whereas the right-hand plot shows the results of the 
best hts reported in Table ^. 



bounds for the branching ratio obtained at 90% confidence level are well above the 
Standard Model prediction of (3.29 ± 0.33) x 10"^. Combining the results obtained for 
the three choices of rrib, we get 

BiB X,7)^c'^EO " = (2.66 ± 0.56e.p t°:l^,J x 10"^ , (33) 

corresponding to a ratio R = Bexp/Bth = 0.81 ± 0.17(exp)loj5(th), which deviates from 
unity by less than one standard deviation. In quoting the theoretical error we assume 
that the spread of the results shown in the three columns of Table ^ represents a rea- 
sonable estimate of the theoretical uncertainty arising from the modeling of the shape 
function. The individual R ratios for the different choices of m;, are displayed in the last 
column in the table. The result (^) is in good agreement with the value (p6|) obtained 
from the extrapolation of the partially integrated branching ratio measured by CLEO 
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Table 5: Results for the total B Xs^ branching ratio, and upper hmits at 90% 
confidence level, obtained from a fit to the CLEO data shown in Table The 
quantity R denotes the ratio of the extracted value for the branching ratio to the 
Standard Model prediction. 





mfo = 4.65 GeV = 4.80 GeV nib = 4.95 GeV 


B{B X,7) [10-4] 
90% CL 

R 


3.09 ±0.66 2.66 ±0.56 2.18 ±0.47 
< 4.66 < 4.06 < 3.26 
0.64/3 0.08/3 0.97/3 
0.94 ±0.20 ±0.09 0.81 ±0.17 ±0.08 0.66 ± 0.14 ± 0.07 



with a photon-energy cutoff at 2.2 GeV. Since the fit to the photon spectrum uses more 
experimental information, we tend to consider (^) to be the more conservative result of 
the two. We will therefore use this value in our further analysis. 

Once high-statistic measurements of the photon spectrum become available, it will 
be possible to determine the parameters of the shape function directly from the data. 
In particular, the average photon energy in B ^ Xg'y decays is a sensitive measure of 
the 6-quark mass. In practice, what can be measured is the average photon energy as a 
function of the cutoff given by 

dB 



Pmin dE^ 

(E,) = — ^ . (34) 

I , „ dB 



/ dE, 



gmin dE^ 



Provided that E"™™ is not too close to the endpoint, this quantity is insensitive to the 
details of the shape function except for the value of mi,. Indeed, at leading order in 
the heavy-quark expansion one simply gets {E^) = mf,/2. The result for the average 
photon energy obtained by including the full next-to-leading order corrections is shown 
in Figure |^. For simplicity, we neglect in this figure the boost between the B rest frame 
and the laboratory frame, which has a very small effect on the average photon energy. 
The different curves in each plot refer to the various sets of shape-function parameters 
considered previously in Figure |^. In the right-hand plot we show the average photon 
energy normalized to mf,/2. We observe that for E"™™ < 1.8 GeV the mean photon energy 
provides a sensitive measure of the mass of the b quark, which to a good approximation 
is independent of other shape-function parameters such as /i^. Asymptotically, for very 
small cutoff values (E^) is lower than mb/2 by about 3%. For i?™" ^ 1.8 GeV, on the 
other hand, the sensitivity to the modeling of Fermi motion quickly increases. For a 
cutoff at 2.2 GeV as employed in the CLEO analysis, there is very little sensitivity to 
the value of rrib. 
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Figure 6: Theoretical predictions for the average photon energy as a function of the 
cutoff -E™™, for different parameters of the shape function. The sohd hnes refer to: 
(1) rub = 4.95 GeV and fil = 0.14 GeV^; (2) m?, = 4.8 GeV and nl = 0.3 GeV^; (3) 
rrib = 4.65 GeV and /u^ = 0.52 GeV^. The dashed hnes show the dependence on the 
value of fi'l for the case = 4.8 GeV, where fi^ = 0.3 GeV^ (solid), 0.45 GeV^ (long- 
dashed), 0.15 GeV^ (short-dashed). The gray lines refer to the Gaussian ansatz for the 
shape function. 



Moments of the type shown in have been considered previously by Kapustin and 
Ligeti employing the heavy-quark expansion but not including the effects of Fermi 



motion. These authors make predictions for the average photon energy and for the width 
of the photon spectrum for cutoff values between 1.8 and 2.2 GeV. Our findings show 
that for i^;™™ > 1.8 GeV there are significant corrections to these predictions caused by 
Fermi motion. This has also been noted by Bauer [^. 

Let us assume that in the future it will be possible to measure the average photon 
energy in B ^ Xg'j decays with a cutoff low enough to be insensitive to Fermi motion 
effects. It is instructive to understand the precise meaning of the 5-quark mass that can 
then be extract from the data. In the limit where Fermi motion can be neglected, it 
follows from (^) that 

1 

/ d?/p?/pPp(?/p) 

(i^.) = f^x^ , (35) 

/ dypPp(yp) 

l-<5p 

where mf, is the pole mass, PpiVp) is the photon spectrum in the parton model, and 
l-6p = 2Ef''/mb. Using that 



1 

/ 



dypVp PpiVp) = Jdz J dyp Ppiyp) + (1 - 5p) j dyp Ppivp) , (36) 
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Figure 7: Theoretical prediction for the function D{5p 



we find for tlie average plioton energy 



(37) 



where 



CE[(ys] 



l-|^ + O(aD-0.97. 

54 TT 



Re[Ct\m,)Cf*{m,)] 



i,j=2,7,8 



(38) 



Tlie prime on tlie sum indicates tliat 7^ (7, 7). The functions dij{Sp) and dyriSp) are 
coUected in Appendix C. The normahzation of ^77 is such that ^77(1) = 0, i.e. apart from 
corrections due to operator mixing the average photon energy without any cutoff is given 
by ^nihCEiois]- Finally, the quantity (5ht in ( P^ ) parametrizes "higher- twist corrections", 
which are neglected in our leading-twist approximation to the shape function. An explicit 
calculation of these corrections gives [pi E3] 



Ai + 3A2 
2m? 



;-0.1±0.4)% 



(39) 



where we have assumed that — Ai = /i^ = (0.3 ± 0.2) GeV^. The contribution of 5ht in 
(^) has a negligible effect. 

In Figure |^, the function D{6p) with 6p = 1 — 2i?™™/mb is shown over a wide range of 
cutoff values. For all realistic values of i?™™ this function takes very small values. The 
corresponding contribution to {E^) in ( P7D is typically less than 1% and thus negligible. 
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Therefore, to an excellent approximation the average photon energy in 5 — * Xs'~) decays 
measures the product of the 6-quark pole mass times the perturbative series CE[oisijnb)]i 
which we have computed to 0{as)'- 

{E,) ^'^CE[as{m,)]^\m\^\ (40) 

It is well-known that the pole mass of a heavy quark is an infrared-sensitive quantity, 
which cannot be unambiguously defined beyond perturbation theory Formally, 



this property appears as a factorial divergence of the expansion coefficients (i.e. an in- 
frared renormalon) in any perturbative series that relates the pole mass to a short- 
distance mass, such as the running mass fnb{mb) in the MS subtraction scheme. On the 
other hand, it is clear that the average photon energy in radiative 5-meson decays is 
a physical observable and as such does not suffer from any ambiguity. Therefore, the 
quantity ml^^ defined in ( |iOD has a short- distance nature. Indeed, it can be shown that 
the renormalon in the pole mass is exactly cancelled by a renormalon (i.e. a factorial 
divergence) in the perturbative series C£;[as] In other words, the short-distance 



mass ml can be related to any other short-distance mass without encountering large 
perturbative coefficients. In particular, we note that 

ml = mb{mb) ( 1 + ^ — ^ • • • 1 • (41) 

Thus, in principle an accurate measurement of the photon spectrum in i? ^ Xg'y decays 
would provide for a theoretically clean determination of the 6-quark mass. In practice, 
this determination will probably be limited by experiment and may not be competitive 
with precision determinations of mb from the analysis of the T spectrum. 



5 Hadronic mass distribution 

In B ^ decays, the invariant mass of the hadronic final state is related with 
the photon energy in the B rest frame through Mfj = m\ — 2mBE^. Therefore, our 
theoretical results for the photon spectrum can be translated into predictions for the 
hadronic mass spectrum. Since experimentally the measurements of the photon energy 
and hadronic mass spectra are quite different, it may be useful to discuss our results also 
in terms of the variable Mh- In the left-hand plot in Figure |, we show the corresponding 
spectra obtained using our standard choices for the parameters of the shape function. 

At this point, it may be worthwhile to recall that the theoretical predictions for the 
photon energy and hadronic mass spectra must be understood in the sense of quark- 
hadron duality. In particular, the true hadronic mass spectrum in the low-mass region 
may have resonance structures due to low-lying kaon states, and thus may look rather 
different from our theoretical predictions. To discuss in more detail how quark-hadron 
duality works in the present case we distinguish two kinematic regions: the "endpoint 
region" and the "resonance region" . The endpoint region of the photon energy spectrum 
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Mh [GeV] Mh [GeV] 



Figure 8: Theoretical predictions for the invariant hadronic mass spectrum for different 
parameters of the shape function 



is characterized by the condition that E^^^ — = 0(A), where A = — mh. It is 
in this region that the effects of Fermi motion are relevant and determine the shape of 
the spectrum. In the endpoint region, the invariant mass of the hadronic final state is of 
order m^A ^ Aqq^, implying that a large number of final states are kinematically ac- 
cessible. Under such circumstances, local quark-hadron duality ensures that the photon 
and hadronic mass spectra are similar to the corresponding inclusive spectra predicted 
by the heavy-quark expansion even without applying a smearing procedure. In the res- 
onance region, on the other hand, the invariant mass of the hadronic final state is of 
order Aqqj-,, implying that the photon energy is very close to the kinematic endpoint: 
^max _ _ O^AQQj-j/m^). The heavy-quark expansion does not allow us to make 
model-independent predictions for the structure of the individual resonance contribu- 
tions. Global quark-hadron duahty can, however, be restored by averaging the spectra 
over a sufficiently wide energy interval, whose size is determined by the average level 
spacing between the resonance states We will see below that in the present case the 
smearing should be done over an interval AM'jj 2GeV^, corresponding to an energy 
interval Ai?^ k, 0.2 GeV. Note that in the case of the CLEO data such an averaging is 
automatically provided by the Doppler shift of the spectrum due to the motion of the B 
mesons produced at the T(4s) resonance, and thus the photon spectrum is expected to 
be dual to the theoretical spectrum over the entire energy range. 

To make these statements more precise, consider the properties of the lowest-lying 
kaon states contributing to S — >■ X<j7 decays, which are collected in Table ^. There 
are six resonances plus a continuum contribution feeding the photon spectrum in the 
energy interval between 2.4 and 2.6 GeV. Hence, an average over this interval should be 
calculable using global quark-hadron duality, although a much finer resolution cannot be 
obtained. In the hadronic mass spectrum, the i^r*(892) peak is clearly separated from the 
rest; however, the next resonances already have widths exceeding the level spacing and 
hence are overlapping. Therefore, we expect that local duality allows us to predict the 
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Table 6: Mean masses and widths of the lowest-lying hadronic states accessible 
in B Xgj decays J5^, and the corresponding photon energies (errors refer to 
changing Mh hy ) 



State H 


Mh [GeVJ 


Vh [MeVJ 


[GeVJ 


K (nvr) 


> 0.629 


continuum 


< 2.60 


ir*(892) 


0.894 


50 


2.56 ±0.01 


iri(1270) 


1.273 


90 


2.49 ± 0.02 


is:i(1400) 


1.402 


174 


2.45 ±0.05 


ir*(1410) 


1.412 


227 


2.45 ±0.06 


/s:*(1430) 


1.428 


103 


2.45 ±0.03 


/r2(1580) 


1.580 


110 


2.40 ±0.03 


fs:i(1650) 


1.650 


150 


2.38 ±0.05 


fs:*(1680) 


1.714 


323 


2.36 ±0.10 


K2{m^) 


1.773 


186 


2.34 ±0.06 



hadronic mass spectrum in the region Mh ^ 1.5 GeV. Indeed, the pattern of resonances 
exhibited in Table ^ suggests a simple but realistic model for the hadronic mass spectrum 
consisting of a single Breit-Wigner peak for the i^'*(892) followed by a continuum above 
a threshold Mcont, which is dual to the higher resonance contributions and given by the 
inclusive spectrum calculated using the heavy-quark expansion. This gives 

dB 2MhNk*'B{B ^ K*^) .dBinci 

± - Mcont) —rj— , (42) 



dMn [Ml - ml, f + ml.T\, ' ^ dM^ 
where 

is the normalization of the Breit-Wigner distribution. The exclusive branching ratio 
for the decay B —>■ K*'~f can either be taken from an independent measurement or 
determined from a fit to the spectrum itself. The continuum threshold Mcont is then 
fixed by the requirement that the total branching ratio be the same as that predicted by 
the heavy-quark expansion, yielding the condition 



Meant , „ 

dM^ TI^ = B(i? X.7)|^ = B{B K*j) , (44) 



/ 



dM, 



H 



where -Econt = 2("^B~^cont)/"^B- order to reduce systematic errors, it will in practice 



be advantageous to normalize both sides of (|4^ ) to the total B Xg^y branching ratio 
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To illustrate this method we consider the central value rrib = 4.8 GeV and take from 
experiment the exclusive branching ratio B{B K*'y) = (0.45 ± 0.17) x 10^'^ [|5ll 



normalized to our fit result for the total branching ratio given in ([33|). This yields 
B{B K*-f)/B{B X,7) = 0.17 ± 0.08. The value of the continuum threshold which 
reproduces the central value of this ratio is Mcont ~ 1.15 GeV, which is close to the 
position of the second resonance i^'i(1270). The dashed line in the right-hand plot in 
Figure ^ shows the hadronic mass spectrum obtained from (^) using these parameters. 
To smoothen out the sharp structures, we take a convolution of this curve with a Gaussian 
smearing function of width ctm^ = 100 MeV, which resembles the binning and resolution 
of a realistic experiment. The result is shown by the solid curve, which exhibits a two- 
peak structure: a narrow peak located at the K* mass, whose width is determined by 
the mass resolution of the experiment, is followed by a broad bump containing a large 
number of overlapping resonances, whose sum is dual to the inclusive spectrum predicted 
by the heavy-quark expansion. The position of the second peak is determined by the 
6-quark mass through M^^"^ ^ (msA)^/^, which is located at 1.6 GeV in the present 
case. This relation is just a reflection of the fact that the 6-quark mass determines the 
average photon energy in i? Xg'y decays (see Section ^). We note that such a two- 
peak structure is indeed seen in the CLEO data on the hadronic mass distribution |l|]; 
however, in view of the large experimental uncertainties it is premature to perform a 
detailed comparison with the data. 



6 Impact of New Physics and constraints on exten- 
sions of the Standard Model 

Measurements of the B —>■ Xg'-f branching ratio have been used extensively to put con- 
straints on the parameters of various extensions of the Standard Model, such as multi- 
Higgs models, supersymmetry, or left-right symmetric models (see Refs. [pO| |3T[| for 
recent reviews). In many of these extensions, the Wilson coefficients of the dipole oper- 
ators O7 and Os in the effective weak Hamiltonian receive additional contributions from 
new flavour physics at a high energy scale. The CLEO measurement of the B Xs'~f 
branching ratio has been used to extract the magnitudes of the Wilson coefficients 
Cj^mw) and Cs{mw) and to compare the results with various model predictions (see, 
e.g., Refs. |5^, |5^). Likewise, new dipole operators of non-standard chirality may be 
induced. In general, such New Physics contributions would enter our analysis through 
non-standard values of the parameters defined in (0) and ([T9|). If these parame- 
ters carry non-trivial new weak phases, there is potential for observing a large direct CP 
asymmetry in the decays B Xg'j, which would be a striking signature for New Physics 



4J]. Here, we shall discuss the impact of New Physics on the CP-average B — > Xs'y 
branching ratio. 

The theoretical predictions for the branching ratio and photon spectrum depend on 
five real combinations of the parameters , as shown in ( ppj ) and (^Tj) . It is convenient 
to introduce the ratio of New Physics contributions to the chromo-magnetic and magnetic 
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dipole operators as 

X = |^. (48) 

A specific New Physics scenario will make a prediction for this quantity. In models 
with dipole operators of non-standard chirality, we assume that (.g/^j is given by the 
same ratio x- Moreover, to simplify the discussion we shall assume that x is a real 
parameter. This is a good approximation whenever there is a single dominant New 
Physics contribution, such as the virtual exchange of a new heavy particle, contributing to 
both the magnetic and the chromo- magnetic dipole operators [Q. With this assumption, 
there are only two independent structures appearing in and (pi]), which can be 
taken to be Re(,^7) and |,^7p + Note that the imaginary part of ^7 enters only in 

combination with the right-handed coupling |^^|, implying that by measuring the total 
branching ratio alone one will not be able to put constraints on the weak phase of ^7. 

As we saw in Section ^, New Physics contributions are unlikely to alter the spectral 
shape of the photon spectrum in the experimentally accessible region. We therefore focus 
our discussion on the total branching ratio and define the ratio 

which directly measures the deviation from the Standard Model. From our result (^) 
it follows that at the level of two standard deviations 0.37 < < 1.25. It will be 
convenient to define a similar ratio of branching ratios for the rare hadronic decays 
B Xsg, which are induced by the chromo- magnetic dipole operator, so that 

^ B{B X,g) ^ |C8K)p + |cfK)p 

Whereas the Standard Model predicts the very small value Bsm(-B Xsg) ~ 0.2%, a 
much larger branching ratio is attainable in models with enhanced b ^ sg transitions 
01^01 • This would increase the production of charmless final states in hadronic B 
decays, which could help to explain the low experimental values of the semileptonic 
branching ratio and charm yield. Although the systematic errors in the measurements of 
these quantities are large, the results favour values of B(_B — > Xgg) of order 10% ^ . 



On the other hand, the CLEO Collaboration has recently presented a preliminary upper 
limit on B{B X^g) of 6.8% (90% CL) The limit is increased to 9.0% if more 



recent charmed baryon and charmonium yields are used ||58[. In our graphical analysis 
below we will present 5% and 10% contours for this branching ratio. Which of these two 
is considered to be a more realistic upper bound is left to the taste of the reader. 
The theoretical predictions for the two ratios R can be written as 

R, = l + A^ix) [M^j) - 1] + A^ix) (l^rr + \^^\' - l) , 

Rg = l + Asix) [M^j) - 1] + A,{x) (l^r + \^7\' - 1) , (48) 
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Table 7: Approximate ranges of x values for various New Physics contributions to 
C7 and Cg, characterized by the particles in penguin diagrams 



New Physics penguins 


X 


gluino-squark 


2-20 


techniscalar 


9 


charged Higgs-top 


1.5-2.4 


Higgsino-stop 


0.2-2.2 


left-right W-top 


0.9 


scalar diquark-top 


-(0.8-1.3) 


neutral scalar-vectorlike quark 


-8 



where 



AM = + XB,. + (1 - x)Ih. ^ ^^^^ ^ ^ ^^^^ 

-D22 + -D27 + -D77 + -D28 + -D78 + -D88 

A,ix) = ^ j^77 + Xi?78 + X^i?88 ^ ^ Q Q25x + 0.0013x^ , 

-D22 + -D27 + -D77 + -D28 + -D78 + -D88 



Mx) = ^^^^ ^ + ^) ^ o.43x(l - X) + 0.50X . 
,1 + r)"^ 



Mx) = jY^,^0.21x\ (49) 

with r = J2i=ihiri""-~^'^^'^^ /C^\'mw) ~ 1.16. The expressions for the functions Ai and 
A2 follow directly from (|20|) , whereas those for ^3 and A4 follow from the result for 
C^\lJ,b) given in (P). Since these functions depend on ratios of coefficient functions, 
their numerical values are rather stable against variations of the input parameters. The 
numerical results quoted above are obtained by taking central values of all parameters, 
and in the case of Ai and A2 using i^;™™ = 1.95 GeV for the cutoff on the photon energy, 
corresponding to the energy interval covered by the data shown in Table § and used to 
derive the result (p3D. 

Approximate ranges of x ior several illustrative New Physics scenarios with 6 — >■ S7 
and b sg penguin diagrams containing new particles in the loop have been discussed 
in Ref. and are collected in Table |^. Our aim here is not to carry out a detailed 



study of each model, but to give the reader an idea of the sizable variation that is 
possible in x- In supersymmetric theories, e.g., penguin diagrams with gluino-squark 
loops imply a positive value of Xj which can be tuned over a large range by adjusting 



^In the notation of this reference x — Csm/^ with ,^gM = —3Cj^{m.w)/Cf^{m.w) ~ —6.05. 
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Figure 9: Regions in parameter space allowed by the experimental constraint on the 
B Xgj branching ratio (bands), and contours for the B Xgg branching ratio 
(dashed lines), assuming Im{^f) = Ct^ = 0. The thin lines in the right-hand plot show 
contours of constant x values as indicated, the black circle corresponds to the Standard 
Model. 



the mass ratio rrig/mg. A detailed analysis of the B X^'-f branching ratio in this 



scenario is presented in Ref. |^6[. Another example with large positive x is provided 
by models with techniscalars of charge | ||5^, |6^, |6l]], whereas an example with large 
negative x is provided by penguin diagrams with new neutral scalars and vector-like 
quarks of charge — |. In general, models with large |x| offer the possibility of having 
a direct CP asymmetry in B ^ Xg'j decays of as much as 10-50% in a large region 
of model parameter space Q]. Examples of scenarios which lead to more moderate x 
values include left-right symmetric models and multi-Higgs models. 

Let us first assume that the New Physics contributions to Cy and Cs do not contain 
new weak phases, and that there are no new operators with non-standard chirality, i.e. 
Im(^7) = C,^ = 0. In the left-hand plot in Figure |^, we show as a function of x the 
allowed ranges for ^7 which satisfy the condition that 0.37 < < 1.25. We also show 
as dashed lines the 5% and 10% contours for the B —>■ X^g branching ratio. There are 
three different regions to distinguish: (a) for x > 3 only positive values Re(^7) < 1 
are allowed, which are close to the Standard Model value ^7 = 1; (b) for — 1 < x < 3 a 
second branch of large negative values of Re (^7) is allowed, which have magnitude several 
times larger than in the Standard Model; (c) for x < — 1 only the first branch remains, 
and the combined constraints from the B Xg'y and B Xsg branching ratios imply 
that — 1 < ^7 < 2.5. The right-hand plot in Figure ^ shows the same information in a 
different way, namely in the plane spanned by the (real) New Physics contributions to 
the Wilson coefficients CY(mw) and Cs{mw)- This figure generalizes the corresponding 
leading-order results discussed in Refs. It is evident that for a given New 

Physics scenario, i.e. for a chosen value of the constraints imposed on the Wilson 
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Figure 10: Regions in parameter space allowed hy the experimental constraint on the 
B — > Xg'^ branching ratio (bands), and contours for the B — > Xgg branching ratio 
(dashed lines), in the general case with non-trivial weak phases and/or non-standard 
chirality operators 



coefficients are quite non-trivial. An example is provided by the "minimal supergravity 
model" 1^ investigated by Hewett and Wells |Q, who perform a scan in the SUSY 



parameter space finding that —2.5 < ^7 < 5.5, and x ~ 1 for choices of parameters 
yielding sizable New Physics contributions to C7 and Cg. From Figure ^ it follows that 
in this model the constraint from the B Xg'-f branching ratio implies —0.1 < ^7 < 1.2, 
which excludes a considerable fraction of the SUSY parameter space, whereas there is 
no constraint from the experimental bound on the B —>■ Xgg branching ratio. Hence, in 
the minimal supergravity model the Wilson coefficients must take values close to those 
predicted by the Standard Model or somewhat smaller. 

The situation becomes more complicated if one allows for complex values of the 
Wilson coefficients arising from new weak phases (i.e. Im(^7) 7^ 0), or considers the 
possibility of having dipole operators with non-standard chirality (i.e. ^7^ 7^ 0) [Q. 
Three illustrative cases are shown in Figure ITO. Note that for large |x| the constraint 



from the B Xsg branching ratio is quite non-trivial and puts a stringent upper bound 
on the combination [Im(^7)]^ -|- l^'f p. 



7 Conclusions 

The inclusive radiative decays B X^'j play a key role in testing the Standard Model 
and probing the structure of possible New Physics. We have presented a detailed study of 
these decays using the operator product expansion for inclusive decays of heavy hadrons 
combined with the twist expansion for the description of decay distributions near phase- 
space boundaries. We have updated (and corrected) the existing next-to-leading order 
analyses of the total B —>■ Xsj branching ratio and added several improvements concern- 
ing the treatment of QED corrections, the analysis of the renormalization-scale depen- 
dence, and the discussion of the sensitivity to New Physics. In particular, we have argued 
that the truncation error of the perturbative expansion in has been underestimated 
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by previous authors by at least a factor of 2. 

Our main focus, however, was to implement a consistent treatment of bound-state 
effects related to the soft interactions of the 6-quark inside the B meson. These effects 
cause the Fermi motion of the heavy quark, which is responsible for the characteristic 
shape of the photon energy spectrum in B ^ Xsj decays. They lead to the main the- 
oretical uncertainty in the calculation of the inclusive branching ratio if a restriction 
to the high-energy part of the photon spectrum is imposed, as is necessary in realistic 
experiments. Fermi motion is naturally incorporated in the heavy-quark expansion by 
resumming an infinite set of leading-twist operators into a non-perturbative shape func- 
tion. The main theoretical uncertainty in this description lies in the value of the 6-quark 
mass. Other features associated with the detailed functional form of the shape function 
play only a minor role, particularly if a partial integration over the decay distributions 
is implied. We have explained how the value of mf, could, in principle, be extracted from 
a measurement of the average photon energy in i? — > Xg'y decays. For the Standard 
Model, we have obtained B{B X,-f) = (2.57 ± 0.26toil) x 10~^ for the integral over 
the high-energy part of the photon spectrum with E}^^ > 2.2 GeV, where the first error 
reflects the uncertainty in the input parameters, and the second one the uncertainty in 
the calculation of Fermi motion. This prediction agrees with the CLEO measurement of 
the same quantity within one standard deviation. From a reanalysis of the CLEO data, 
we have obtained values for the total branching ratio that are consistent with the Stan- 
dard Model prediction of (3.29 ± 0.33) x 10~^. In the future, an effort should be made 
to lower the cutoff on the photon energy to a value of 2 GeV or less, even if this would 
increase the experimental systematic errors. The benefit of such a low cutoff would be 
that the calculation of the branching ratio becomes insensitive to the effects of Fermi 
motion, reducing the theoretical uncertainty to the level of 10%. 

Besides the photon spectrum, we have studied the invariant hadronic mass distribu- 
tion in radiative B decays. Investigating the pattern of individual hadron resonances 
contributing to the spectrum, we have argued that a complementarity between the in- 
clusive theoretical distribution and the true spectrum should set in for Mh ^1-5 GeV. 
This leads us to a simple description of the hadronic mass spectrum with only a single 
parameter, the B — > K*j branching ratio, to be determined from experiment. 

Finally, we have investigated the sensitivity of the B — >■ branching ratio and 
photon spectrum to New Physics beyond the Standard Model and set up a formalism 
to include possible non-standard contributions in a straightforward way. New Physics 
contributions enter our predictions through the values of parameters which are 
defined in terms of Wilson coefficient functions at the scale mw- This formalism allows us 
to account for non-standard contributions to the magnetic and chromo-magnetic dipole 
operators, as well as operators with right-handed light-quark fields. In the context of 
a particular model, all that is needed is to perform a matching caknilation determining 
the Wilson coefficients in the effective Hamiltonian at the weak scale. We find that, 
quite generally. New Physics contributions would not affect the shape of the photon 
spectrum, but of course could change the total branching ratio by a considerable amount. 
This implies that the analysis of the photon energy and hadronic mass spectra, which is 
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crucial for the experimental determination of the total branching ratio, can be performed 
without assuming the correctness of the Standard Model. On the other hand, the total 
branching ratio will provide a powerful constraint on the structure of New Physics beyond 
the Standard Model, as we have illustrated with some specific examples. 

We believe that our work eliminates the remaining elements of model dependence 
present in previous studies of the photon spectrum in B ^ Xg'j decays and, therefore, 
provides a firm theoretical basis for analyses of experimental data on inclusive radiative 
B decays. We are confident that in the near future, when measurements of the photon 
spectrum with high statistics will be performed, it will be possible to derive a value 
for the B Xg'j branching ratio that is less model-dependent than existing ones, thus 
providing one of the most sensitive tests of the flavour sector of the Standard Model. 

Note added: While this paper was in preparation the CLEO Collaboration presented 
a preliminary update of the B —>■ X^j branching ratio, yielding B{B —>■ Xg'-f) = (2.50 ± 
0.47 ±0.39) X 10"^ H. This value has been obtained using the original analysis adopted 



in Ref. [|^. We expect that using our improved theoretical predictions the central value 
will increase to about 2.8 x 10~^, which is close to the prediction of the Standard Model. 
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Appendix A: QED corrections to Cr{mij) 

QED corrections affect the theoretical prediction for the B Xs'~f branching ratio in 
three ways: there are 0(a) matching corrections to the Wilson coefficients Ci{mw) at the 
weak scale, there are 0{q) contributions to the matrix elements of the operators in the 
effective Hamiltonian at the scale /i;,, and there are 0{aL) corrections to the evolution of 
the operators from the scale mw down to the scale /i;,, where L = ln{mw / fih) . The latter 
corrections are logarithmically enhanced. They can be accounted for by including the 
0{a) contributions to the anomalous dimension matrix of the operators in the effective 
Hamiltonian in the solution of the renormalizat ion-group equation. On the other hand, a 
complete calculation of the remaining 0{a) corrections would be extremely cumbersome. 
Fortunately, it is likely that these corrections will be smaller than the next-to-next-to- 
leading QCD corrections of order a^. 

When QED corrections are included, as many as twelve operators in the effective 
Hamiltonian mix under renormalization. Besides the current-current operators Oi and 
O2 and the dipole operators O7 and Og, these are four QCD penguin and four electroweak 
penguin operators. However, it turns out that to a very good approximation the mixing 
of (Oi, O2, O7, Og) with the penguin operators can be neglected. This approximation 



34 



reproduces the leading-order results in (||) exactly but for the terms J2i hif]"-^ and hii]"-^ 
in Cj'^ and Cg'\ Numerically, Cj^\mi,) is reproduced with an accuracy of 3 x 10~^, and 
Cg'\'mh) with an accuracy of 3.7%. 

Including QED corrections, the one-loop anomalous dimension matrix in the trun- 
cated operator basis is 

(A.l) 



7 = — 7o + — To , 



47r 



47r 



where 



7o 



-2 
6 



6 3\ 

O 416 70 
81 27 



f 

V o-f fy 



^-1 



To 



V 



208 


8 
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81 


27 


208 


116 




243 


81 




16 


8 




9 


3 







8 
9 


/ 



(A.2) 



We have determined the entries of the QED matrix To using the results of Ciuchini et 
al. for the contributions to each entry of the QCD matrix 70 taking into account 
their color structure and the electric charges of the quark fields involved. We collect the 
Wilson coefficient functions into a vector C and write their evolution as 



C{[i) = U(/i, mw) Cimw) , 
where the evolution matrix satisfies the renormalization-group equation 



(A.3) 



-T- U(/i, mw) = 7^ U(;U, mw) ■ 
djj 

To solve this equation, we make the ansatz 

U(/i, mw) = K(/i) Uo(/i, mw) K'^{mw) 



(A.4) 



(A.5) 



where Uo(/i, mvi/) is the leading-order QCD evolution matrix. To first-order in the 
electro-magnetic coupling a we may write 



K(/i) = l + a 



Kq 

asifi) 



(A.6) 



where the ellipses represent terms that do not contribute at leading-logarithmic order. 
Inserting this ansatz into the renormalization-group equation satisfied by K(yu), we ob- 
tain Kq from the solution of the algebraic equation 



2/3oKo = rJ + [7(^',Ko], 



(A.7) 



where /3o = y is the first coefficient of the QCD /3-function for Uf = 5 light quark fiavours. 
Let us introduce the matrix V that diagonalizes the QCD anomalous dimension matrix. 



I.e. 



V So V = diag(7i, 72, 73, 74) 
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and denote k = V ^Kq V and g = V ^Fq V. Then the solution of ( |A.7| ) yields 



^^^■"2/3o-(7.-7,)- ^^-^^ 

The result for the evolution matrix U(/i, mw) can now be written in the form U(/i, mw) = 
VuV^-*^, where 

This general result has previously been derived by Buchalla et al. [0. The evaluation 
of this relation for our particular case leads to the expression for the QED coefficient 
C'f™\lJi,h) given in ([TT|). The coefficients of C'f'\mw) and cf^{'mw) in this result are 
exact, whereas the remaining terms are only approximate because of the truncation of 
the operator basis. However, as in the case of the QCD evolution we expect that from 
a numerical point of view the truncation of the basis is justified. Expanding the result 
(pAp to first order in as we recover the formula of Czarnecki and Marciano ||33|1 : 



— ctH^.) = f - ^ Crim^)] In ^ + . . . . (A.ll) 



47r V243 9 ' ' V /if, 



Numerically, the resummed expression in (jl^) is smaller than the naive result ( |/rTT| ) by 
a factor of about 0.55. 

We note that Czarnecki and Marciano also include a particular type of matching 
correction to the coefficient Cj{mw) at the weak scale, arising from fermion-loop in- 
sertions on the ly-propagator in the top-quark-VT penguin diagram [Q. The result- 
ing contribution is AC7(mvi/) ~ 0.53q;, which at the scale rri}, yields a contribution 
/S.C-j{mij) 0.36a ~ 2.6 x 10~^. Since there are many other matching corrections that 
have not yet been calculated, and since there are similar 0{a) contributions to the ma- 
trix elements of the local operators Oi that are neglected, we see no compelling reason 
to include this particular matching contribution. Its effect on the total branching ratio 
does not exceed the level of 1% and is thus safely within the theoretical uncertainty of 
±2%, which we assign to higher-order electroweak corrections. 



Appendix B: Photon spectrum in decays of B mesons 
produced at the T(4s) resonance 

We denote by E^^"^ the photon energy measured in the laboratory, and by the one 
measured in the B rest frame. If dB/di?^ is the photon spectrum in the B rest frame, 
the corresponding spectrum measured in the laboratory is 

dB 1 /• ,„ 1 dB 



dM^b (3+ -(3- J ' E^ dE^ 



d^7TrTTr, (B.l) 
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where 



\ah\ 
7 ' 



(B.2) 



The maximum photon energy in the laboratory is {E, 



\ab\ 



/3+E 



max 

7 



It is straightforward to calculate from ( |B.1| ) the effect of the boost on the partially 
integrated branching ratios. Let us define the difference 



A(^o) 



y-^-y jmax. ^-y 

7 d^lab 



dE^ 



Eq ' Eq 

Then, provided that Eq < we find that 



dB 



(B.3) 



A{Eo) 



P+Eo , . Eo , 

1-P y dB / /?+Eo\ 1 + / dB / 



Eo 



P-Eo 



E^ 



(3'^El (d_^_m\ 

6~ \dE^dE) 



(B.4) 



E=Eo 



i.e. the effect is quadratic in the small quantity (3. With = 2.2 GeV as in the CLEO 
analysis, we find that A(i?o) ~ —(0.04-0.09) x 10""^ depending on the parameters of the 
shape function, which is a very small effect. 



Appendix C: Results for the functions Sij{y) and dij{S) 

The functions Sij{y) entering the theoretical expressions for the photon spectrum are 
given by the derivatives of the functions fij{S) in ([T6|) through Sij(y) = f[j{l — y). 
Explicitly, we find 



sn{y) 

S78{y) 



1 r 



7 + y-2y^-2il + y) ln(l - y) 



1 / 2(2-2y + y^ 
27 1 y 



ln(l - y) + 21n 



y^ 



-2y'-y 



(1-1/) 



9 



16 r 

522(1/) = ^ J dx{l-x 


S27(y) = -Ss28{y) = I dxRe 



ln(l -y) + l + ^ 

y 4 



X \zj 2 



(C.l) 



37 



The functions dij{S) entering the expression for the average photon energy in (p7|) 
are given by 



dij{S) = dx fij{x) - 5 fij{5) . 



We find 

^77(5) 
C?88(5) 



^78 (5) 



4^2 ^4 

(3-5n(5 \ 

9 ^ ' 3 27 6 



27 



In 



52 53 

ln(l-5) + 5 + - + - 



+ 



27 
6 



r ^ '^^ 



-L2(l-5)+ (5 + ^)ln5-5 



ln(5-(5 

752 
- -r- + 



53 



4 36 ^ 8 



16 



27 



1 

y"da;(l-a;)(l-X5)^ 



+ 



1 



(C.2) 



d^,{6) = -Sd^siS) = — jdx(l- xsf Re 



G 



+ 



X 

2i, 



The function ^77 (5) is defined as 



rf77(5)=rf77(5)+Ml + ^ln5)+^. 



(C.3) 



(C.4) 



The second term is the contribution of the Sudakov logarithms. It is this contributions 
which, in large orders, develops a factorial divergence that cancels the infrared renor- 
malon of the pole mass The last term in ( |C.4| ) is adjusted such that dTj{l) = 0. 
This condition defines the coefficient C^fas] in (|37|). 
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